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We study theoretically the phase diagram of the double- 
layered perovskite manganites taking into account the orbital 
degeneracy, the strong Coulombic repulsion, and the cou- 
pling with the lattice deformation. Observed spin structural 
changes as the increased doping are explained in terms of the 
orbital ordering and the bond-length dependence of the hop- 
ping integral along c-axis. Temperature dependence of the 
neutron diffraction peak corresponding to the canting struc- 
ture is also explained. Comparison with the 3D cubic system 
is made. 
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I. INTRODUCTION 

Perovskite man- 
ganites (lili_ x , A x )„ + iMn n 3n+ i (R=La, Pr, Nd, Sm ; 
A= Ca, Sr, Ba ; n — 1,2, oo) have recently attracted 
renewed interests from the viewpoint of the .-close con- 
nection between magnetism and transport E a The sta- 
bility of the perovskite structure enables the preparation 
of high quality single crystals with systematically chang- 
ing the carrier density and the bandwidth by control- 
ling R and A atoms. In addition, the crystal structure 
changes from two-dimensional (n — 1, single layer) to 
three-dimensional (n = oo, cubic) as n increases. The 
dimensionality of the electronic structure is controlled 
also by the spin and the orbital orderings, for example, 
the transfer along c-axis is forbidden by the anti-parallel 
spin configuration and/or the d x 2_ y 2 orbital orderingE'u 
Therefore even in the isotropic (n = oo, cubic, 
113-system) crystal structure, the electronic dispersion 
can be quasi- two-dimensional llayered antiferromagnetic. 
(spin A) in P^Sr.MnOaJita Ndi_ x Sr x Mn0 3 £ra 
Lai-zSr^MnOsli-j) or quasi-one-jdimensional (rod type 
AF (spin C) in Ndi^Sr^MnOaBB). Another impor- 
tant related issue is the spin canting. As discussed by 
de Gennesc£l, A-type antiferromagnet near the parent in- 
sulator is unstable toward the spin canting when holes 
are doped. This is because the kinetic energy gain along 
c-axis wins the energy cost of the AF exchange interac- 
tion for small canting angle. In the metallic A-type AF 
phase mentioned above, however, no spin canting has 
been observedfl 

Double-layered manganites with n = 2 (327-system) 
offer an interesting opportunity to study the interplay 
of the spin, the orbital and the crystal structure in 
the above-mentioned issues of dimensionality and spin 



canting. With increasing x, the spin ordering within 
one double-layer changes from the ferromaapetic. one 
(spin F, x = 0.3) to the spin A (x = 0.5)M-B"N For 
0.4 < x < 0.48 the diffraction peak observed in the 
neutron-scattering experiments indicates the coexistancc 
of the spin F- and A-components at low tempera turns,, 
which has been interpreted as the spin cantingOt3E§E£l 
However it could be interpreted also in terms of the 
phase separation, and further theoretical studies are 
needed. Lattice, deformation depending upon x is also 
observed,E3E3c3 implying active contribution to the spin 
transition, in comparison with 113-system. In 113- 
system, there is the one to one correspondence between 
the crystal deformation and the orbital ordering, e.g., 
the c-axis contracts in the metallic spin A statecj where 
the d x 2_ y 2 ordering occursJj In 327-system, on the other 
hand, MnOg octahedra elongates along c-axis during the 
spiij, structural change from spin F to cant, and to spin 
AB3 This elongation decreases with increasing x upto 
x = 0.5 continuously, bat even in the spin A phase a lit- 
tle elongation remains,c-3 though this deformation prefers 
the d^ z 2_ r 2 . 

In this paper, we study the spin and the orbital order- 
ing in the double-layered compounds. Possibility of the 
spin canting is studied for general x, which includes the 
de Gennes's canting mechanism as a special case. In the 
metallic region, the spin A is locally stable against the 
spin canting. Spin transition from the spin A is there- 
fore discontinous one to the spin canting or to the spin F, 
depending on the ratio of the transfer integral and the su- 
perexchange interaction between ti g spins. For the spin 
canting to occur in the metalic region, the small trans- 
fer integral along c-axis, i.e., the planer orbital ordering, 
turned out to be indispensable. The canting angle is sen- 
sitive to the bond-length along c-axis. With this mecha- 
nism, calculated mean-field phaaejdiagram together with 




can qualitatively 
and T-dependencajp , 



the observed lattice dis 

explain the observed x-UO and 1 -dependencftp , pi 
the spin ordering and also the spin-anisotropyE3t3o"EII 
in 327-system. 

The plan of this paper follows. In Sec. II, the model 
and the formulation are given. Results and discussion are 
presented in Sec. III. Notations used here are standard 
as in ref. 8. 



II. MODEL AND FORMULATION 



We start with the Hamiltonian 
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where 7 [= a(d x 2_ y 2), 6(<i3 Z 2_ r 2)] specifies the orbital and 
the other notations are standard. The transfer integral 
tjy depends on the pair of orbitals (7, 7') and the direc- 
tion of the bond (i, j)E3, Jh is the Hund coupling between 
e g and t% g spins, and Js is the AF coupling between near- 



est neighboring t 2g spins. H a 



represents the on-site 



Coulomb interactions between e g electrons, and is given 



characterized by an angle 77 which is (tt) for spin F 
(A). As for the orbital degrees of freedom, we consider 
two sublattices I, and 17, on each of which the orbital is 
specified by the angle Qiji asO 
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We also consider four types of orbital-sublattice ordering, 
i.e., F-, A-, C-, G-type in the cubic cell. Henceforth, we 
often use a notation such as spin A, orbital G (0j,0jj) 
etc.. Denoting the wave vector of the spin (orbital) or- 
dering as qs (qr), the ground state energy is given as a 
function of the spin ordering (77, qs), the orbital ordering 
{@iji, Qt), and the lattice distortion (g, r, v). 



III. RESULTS AND DISCUSSIONS 



(2) A. Phase diagram without canting and lattice 

distortion 



where the spin operator for the e g 



electron is defined 



as Sp 



■ I d[ dapd^p with the Pauli matrices <?, 

while the orbital isospin operator is defined as Tj = 
\ £ d\ a^i d^i a . Coefficients of the spin and isospin 
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operators, i.e., a and (3, are given bycrEa a — U — ^ > , 
and P = U — kf > . The minus sign in Eq. (||) means 
that the Coulomb interactions induce both spin and or- 
bital (isospin) moments. The parameters a, [3, to, used 
in the numerical calculation are chosen as to = 0.72 eV, 
U = 6.3 eV, and Jq^ 1.0 eV, being relevant to the actuaJ 
manganese oxides.Qta The electron-lattice interaction isa 



#ci- P h = + \g\r 



■Ti 



(3) 



where g is the coupling constant and r (ili) is the mag- 
nitude (direction) of the lattice distortion of the Mn06- 
octahedra. Values of r and v are taken from the ob- 
served elongation along c-axis in (Lai_ x , ST x ) 3 Mn 2 0i 
(0.3 < x < OAJB leading to r - 0.01, v//z. We calculate 
the ground state energy by the jpeanfield approximation 
with the two order parameters ,□ 



rS = ( S, ... ) + 2| (^2 S 



T 



(4) 
(5) 



As the exchange interaction between two double-layers 
is reported to be less than 1/100 compared with the intra 
double-layer one,E3 we consider an isolated double-layer, 
for which the Brillouin zone contains only two fc-points 
along c-axis. We consider four kinds of the spin alignment 
in the cubic cell: spin F, A, C and G (NaCl-type). For 
spin A, we also consider the possibility of the canting 



FIG. 1. 

Figure, [l] shows the phase diagram at zero temperature 
in the plane of the doping concentration (x) and the AF 
interaction between t 2g spins (Js), without the electron- 
lattice coupling (<? = 0). Orbital shape is represented 
by p- (planer, orbital F (0,0), i.e., d x i_ y i) and n- (non- 
planer, orbital G (100,-100) around x = 0.1 and orbital 
A (-40,-40) around x — 0.7) hereafter. The possibility 
of the spin canting is not taken into account, i.e., rj = 
or tt. For Js = 0, spin F is the most stable for almost 
the whole region of x, except around x = 0. This fer- 
romagnetism comes from the superexchange interaction 
between e g spins for small x [x ~ 0.1), and from A he. 
double-exchange interaction for larger x (x ~ 0.7)J§E!l 
Nonmonotonic behavior with double peaks (x = 0.1, 0.7) 
of the phase boundary between spin F and A (Jg(FA)) 
is then attributed to the cross-over from the super- to 
the double-exchange interaction, similarly to the case of 
113-system.El We identify the spin transition from spin 
F to spin A with inqceasitig x observed experimentally 
for 0.3 < x exp , < 0.5QO~c3 with that around x ~ 0.15 
when Js ^ 0.006 eV (for example, arrow (a) or (b) in 
Fig. 0). This value of J5 is roughly of the same order 
in magnitude as the 113-caseJ3lj Depicted orbital align- 
ment is the one optimized at each point of the phase 
diagram. For most of the phase diagram, d x 2_ y i is the 
most stable one reflecting the two-dimensional nature of 
the double-layered structure. The spin C phase com- 
pletely disappears. One might regard this self-evident 
because of the two-dimensionality. This, however, is not 
so trivial because the valency along c-axis can still lower 
the energy by the bonding/anti-bonding splitting for the 
double-layered structure. 
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B. Spin canting 

We now study the spin canting between spin F and A. 
First we consider the possibility of the spin canting with 
fixed Js = 0.004 eV, shown as the arrow (a) in Fig. [I]. 

FIG. 2. 

The x-dependence of the optimized canting angle 77 is 
shown in Fig. |^. The electron-lattice coupling is not 
considered (g = 0). The canting angle changes continu- 
ously from 7r (spin A, x = 0) to zero (spin F, x ~ 0.1), 
corresponding to the spin canting around the parent in- 
sulator. In this region, the orbital alignment almost re- 
mains that for x — (orbital G (100,-100)), implying 
that the ordering is still dominated by the superexchangc 
interaction in this region. This is consistent with the in- 
terpretation that the peak of the phase boundary Js (FA) 
around x — 0.1 is attributed to the superexchange inter- 
action. Coexisting double-exchange interaction due to 
the introduced carriers competes with Js, leading to the 
spin canting seen for < x < 0.1 in Fig. || through the 
conventional mechanism a la de GennesEJ for 113-system. 
Battle et al. observed the coexistance of the spin F and 
A by the powder diffraction of the poly-crystal sample 
La 2 SrMn 2 7 [x = 0).EI This can be explained by small 
amount of carriers introduced by the oxygen deficiencies 
which brings about the spin canting by this mechanism 
for small x region, as reproduced in Fig. |. When x 
goes beyond x ~ 0.125, 77 discontinuously changes from 
zero to 7r, corresponding to the first-order transition from 
spin F to spin A. This behavior without canting is due 
to the discontinous orbital transition. With increasing 
the double-exchange interaction for x > 0.1, the orbital 
discontinuously changes from orbital G (100,-100) to or- 
bital F (0,0), i.e., d x 2_ y 2 1 in order to maximize the ki- 
netic energy gain. With fixed Js as in Fig. ^, the spin 
transition between spin F and A is accompanied with 
this orbital transition (n-spin F to p-spin A). There- 
fore, the spin canting obs«iV|i 
the single-crystal samples)! 
reproduced in Fig. g. 

In order to explain this experimental observation, we 
look for the possible mechanisms for the continous change 
of spin/orbital structure starting from the p-spin A. One 
possibility is to change the orbital wayerfurLction via the 
coupling to the lattice. In 327-systemE2lE2rE£l the lattice 
constant c is longer than a, b, which prefers d^, z 2_ r 2 com- 
pared with d x 2_ y 2, and c changes as x. This corresponds 
to the "magnetic field" to the orbital pseudo-spin, and 
tends to enhance the transfer along z-direction and hence 
the spin canting. 

FIG. 3. 

In Fig. [3] shown the phase diagram in the plane of x and 
gr (<? : coupling constant, r : elongation along c-axis). 



,0.4 < x exp , < 0.5 in 
6 ' therefore, cannot be 



This gr corresponds to the "orbital magnetic field" to 
prefer d 3z 2_ r 2. Here Js is fixed to be 0.004 eV whictuk 
comparable to the value appropriate for 113-system. Em 
The p-spin A state (shaded in Fig. ||(a)) is surrounded 
by the discontious "orbital spin-flop" transition. This 
is because the "orbital magnetic field" gr is applied not 
perpendicularly but anti-parallel to the c? ;c 2_j / 2-direction. 
Therefore it seems unlikely that the spin canting occurs 
between the p-spin A and spin F with different orbital 
structures, e.g., the arrow (a) in Fig. |l|. 

We now look for the possibility of the spin canting with 
the orbital structure fixed to be d x 2_ y 2 , i.e., the arrow (b) 
in Fig. [I]. The corresponding Js value might appear to 
be too small (~ 0.0002 eV), but Js is very sensitive to 
the bond- length I ( Js ~ l~ 14 as discussed later), and the 
lattice elongation along c-axis suggests smaller value of 
Js in 327-system compared with 113. Furthermore, the 
x-dependence of I induces the x-dependence of Js and 
the arrow (b) in Fig. |l| has finite vertical component. 
Another remark is that even though the orbital pseudo- 
spin Ti is fully polarized along d 2 .2_ 2/ 2-direction, the wave- 
fuction is the hybridized one with d x 2_ y 2 and d 3z 2_ r 2. 

Namely the d 3z 2_ r 2 has the weight ~ (to/ff) , which 

induces the effective transfer t z between layers. However 
t z is much smaller than to, which is crutial to the spin 
canting as discussed below. 

FIG. 4. 

In Fig. H shown the energy as a function of the angle 77 
between the spins of two layers for several values of Js- 
Spin F for Js — is due to the above-mentioned double- 
exchange interaction with the weig ht ~ (to/ffy . For 

Js = 0.42-0.5 meV, the optimized angle 77 is the inter- 
mediate between 77 = and r\ = 7T, and the canting state 
is realized. Therefore we conclude that the spin canting 
occurs between p-spin A and p-spin F states when Js 
changes. 

FIG. 5. 

In order to study further this canting state, we introduce 
a simplified model shown in Fig. ^. We assume that the 
spin polarization is perfect, i.e., all the electron spins are 
aligned along the ordered direction. Therefore we con- 
sider the spinless fermions. The density of states is taken 
to be a constant Np for each layer. We take the hole- 
picture, and x holes are occupying this constant density 
of states from the bottom. These two bands are split 
into bonding and anti-bonding ones with the splitting 



A = t z cos 2 = t z (,S3 The kinetic energy gain AEk in (^) 



is given by 
AE kin (£)< 



-tl-N F -e (far £ < & = j^) 
-t x -x-£ (for £>£ c ) 



(7) 



while the energy cost of the exchange interaction is 
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Jscos r\ = J s (2£ 2 - l) . The lower line of Eq. (0) is ob- 
tained hv de Gennes, and if this holds the canting always 
occursila The new aspect here is that AEki n (£,) oc £ 2 
when the splitting A = t z {; is smaller than the Fermi en- 
ergy e F = x/N F and both the bonding and anti-bonding 
bands are occupied. Therefore the spin A structure 
(£ = 0, rj = 7r) is at least locally stable when 2Js > t 2 N F . 
This condition can be satisfied when the orbital is al- 
most d x 2_ y 2 and t z is much reduced from to- For general 
orbital configuration t 2 N F is order of magnetude larger 
than Js ■ Now we look for the optimized £ to minimize the 
total energy AE(£) = AE km (£) + AE ex (£) . When £ c > 1 
(x > t z N F ), only the upper line of Eq. (Q) is relevant 
and AE = (2J S - t 2 z N F ) ■ £ 2 . Therefore £ jumps from 1 
(spin F) to (spin A) as Js increases across t 2 N F /2. 



FIG. 



When £ c < 1, the optimized £ as a function of Js 
is given in Fig. |^. As J5 increases, the spin structure 
changes as spin F (Js < t z x/4) —> spin canting (t z x/4 < 
Js < t 2 N F /A) — > spin canting with fixed canting angle 
(t 2 z N F /A < J s < t 2 z N F /2) -> spin A (t 2 z N F /2 < J 8 ). 
Note that the canting angle continuously evolves from 
spin F, but jumps at the transition to the spin A. Sum- 
marizing, x < t z N F and t z x/4 < J s < t 2 N F /2 are the 
condition for the occurence of the spin canting. Consider- 
ing N F ~ 1/ioj and Js is about two orders of magnitude 
smaller than to, this condition is satisfied only when t z 
is much reduced from to and x is small. Spin canting 
state does not appear in Fig. || because of the large t z 
for n-spin F. 

Spin ca: 

Xexp. ^ - 51 — 

in this region. This is consistent with the spip_ea 
observed within the a6-plane in this regionEj'LJ 
cause the planer orbital leads to this anisotropy of the 
easy axis according to the spin Hamiltonian derived by 
MatsumotoLJ taking into account the spin-orbit interac- 
tion. The observed spin transition from F to A with 
increasing x should therefore corresponds to that from 
p-spin F to p-spin A in Fig. ^, for which the increase in 
Js with increasing x is needed (arrow (b) in Fig. [l]). 

This x-dependence of Js can be explained in terms 
of the bond-length dependence of the inter-layer hooping 
integral. According to the pseudo-potential theoryEJ, the 
hopping integral between the d orbitals depends on the 

bond length / as t oc l~ 7 , leading to £0 tf 9 / (^2,) <x 

l c , where l c denotes the bond-length along c-axis. The 
observed lattice contraction along c-axis with increasing 
xp3 therefore, gives rise to the decrease of £0, leading 
to the spin structure changes from spin F to cant, and 

(**2g) x mcrea ses due to this 



rved in the metallic region 0.4 < 
6 ' implies therefore the planer orbital 



taxis 
3 be- 



to spin A. J^ oc 

lattice contraction along c-axis, which can explain the IE- 
dependence of Js needed for the change from p-spin F 
to p-spin A, as shown by the arrow (b) in Fig. 



The neutron-diffraction experiments have revealed 
that the samples with the spin canting at zero temper- 
ature (0.4 < frrpp. < 0.48) become the spin A at higher 
temperatures.E3 This can also be explained in terms of 
the lattice deformation ^Jiallpws. The observed lattice 
contraction along c-axisE3tjE3 with increasing tempera- 
ture brings about the decrease of £0 oc t z j Js oc ij,, leading 
to the spin transition toward the spin A discontinously. 
Hirota et al. also attributed thiSplC-dependent spin tran- 
sition to the lattice deformation.Ej 

As seen above, the observed lattice deformation plays 
an important role to explain the spin transition. In our 
scenario, the bond-length dependence of the inter-layer 
hopping is the most important mechanism which relates 
the lattice- and the electron-systems. One might, how- 
ever, attribute it rather to the Jahn- Teller type electron- 
lattice interaction through the electro-static coupling (in 
other word, the "orbital magnetic field"). In this sce- 
nario, the lattice contraction along c-axis reduces the hy- 
bridization of d 3z 2_ r 2 orbital as x increases—namely the 
stabilization of d x 2_ y 2 which prefers spin A.t3 According 
to this scenario, the temperature dependence of the spin 
canting mentioned above is alsp,explained by the change 
in the hybridization of d 3z 2_ r 2t3 which is brought about 
by the temperature-depenedent lattice elongation. This 
scenario, however, results in the decrease of the inter- 
layer hopping as x (T) increases, being opposite to our 
scenario, where little change in the orbital ordering and 
the increase of the inter-layer hopping t z oc l~ 7 during 
the spin transition are predicted. Main distinction be- 
tween these two scenarios is whether such a large change 
in the hybridization of d 3z 2_ r 2 occurs during the spin 
transition or not. As seen in Fig. || and Fig. |3[ however, 
the spin transition should be of the first order without 
canting if it is accompanied with the large change in the 
orbital. Experimentally, the spin easy axis falls ontc 
plane discontinously around x exp = 0.32 (spin FjEj 
and the spin transition from F to cant, and to spin A 
occurs with this easy axis being unchanged. This behav- 
ior implies that the hybridization of d3 Z 2_ r 2 decreases 
discontinously before the spin transition, and during the 
transition the wave- function is almost planer as d x 2_ y 2, 
namely the transition is between p-spin F and p-spin A 
with canting. Though, even in this case, the "orbital 
magnetic field" brings about a little change in the hy- 
bridization of d 3z 2_ r 2 via the change of the weight men- 
tioned before, (to/fij — * (to/ {ft — 5"")) , it is a minor 
correction comparing with the sensitively changing bond- 
length dependence, t z oc l~ 7 , Js oc l^ 14 - 

For the spin canting, £0 = t z x/AJs < 1, and hence the 
planer orbital during the spin transition turned out to 
be essential. In 113-system, the spin F phase neighbor- 
ing the metallic spin A has the quasi-three-dimensional 
orbital ordemig for the realistic parameters in the mean 
field theory.LrB Therefore the canting does not occur be- 
cause the orbital structures are different between spin A 
and F, and the transition between them is always discon- 
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tinous. In the real 113-system no anisotropy has been re- 
ported in the spin F state, and there is a theoretical sug- 
gestioruof the orbital liquid state for the spin F metallic 
region. E3 In any case the orbital state is different from the 
d x 2_ y 2 in the metallic spin A state. In 327-case, the lay- 
ered crystal structure brings about the two-dimensional 
orbital anisotropy, reflecting the interplay of the dimen- 
sionality of the crystal structure and that of the orbital 
ordering. 

FIG. 7. 

Figure |?] summarizes our interpretation of the exper- 
iments (notations Xth. and x exp , are used to prevent 
the confusion due to the quantitative discrepancy of x 
between the calculations and the experiments). The 
hatched region corresponds to the spin canting phase due 
to the planer orbital for the p-spin F. J$ (FA) is the 
phase boundary with the first-order transition between 
spiiLjF (or canting) and spin A. Spin F phase at x exp , — 
0.3E3 corresponds to the n-spin F phase. With increasing 
x, spin F changes from n- to p- with the orbital transition 
into d x 2_ y 2 around xth. ~ 0.125. This transition is of the 
first-order, leading to the discontinous reorientation of 
the spin easy axis from being along c-axis into within the 
a6-plane by considering the spin-orbi t int ccaction.a This 
explains the observed spin anisotropyE3c3'E3'EZl where the 
easy axis, pointing parallel to the c-axis for x exp . = 0.3, — 
discontinuously turns onto the a6-plane at x exp . ~ 0.32,E£ 
for T < 100 K. The p-spin F phase around x th . > 0.125 
cants readily within the a6-plane because of the decrease 
in £o oc l 7 c due to the contraction of the bond-length 
along c-axis with increasing x, and then discontinously 
jumps to the p-spin A. corresponding to the observation 
for 0.4 < x exp . < 0.5.li3iaM 

In conclusions, we have studied the spin and the or- 
bital ordering in 327-system, taking the orbital degrees 
of freedom and the strong Coulombic repulsion into ac- 
count. Observed ^-dependence of the spin ordering and 
its anisotropy were qualitatively explained in terms of 
the orbital ordering and the observed lattice deforma- 
tion along c-axis. The observed temperature depencence 
of the spin structure is also explained by the lattice defor- 
mation. The difference between 327-system with canting 
and 113-system without canting is understood in terms 
of the difference in the dimensionality of the crystal stru- 
cuture. 
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Figure captions 

Fig. [l| Phase diagram at zero temperature in the plane 
of the doping concentration (x) and the AF interaction 
between ti g spins (Js), without electron-lattice coupling 
(g = 0). Orbital shape is represented by p- (planer, or- 
bital F (0,0), i.e., d x 2_ y 2) and n-(not planer, orbital G 
(100,-100) around x = 0.1 and orbital A (-40,-40) around 
x = 0.7). The possibility of the spin canting are not 
taken into account (77 = 0, 7t). Depicted arrow (a) ((b)) 
corresponds to the spin transion with fixed (increasing) 
Js- 

Fig. |^. The ^-dependence of the optimized canting angle 
77 with g = 0. 

Fig. ||. The phase diagram as a function of the dop- 
ing concentration x and the electron-lattice coupling g. 
p-spin A is always surrounded by the discontinous bound- 
ary with respect to r\. 

Fig. |i[ yy-dependence of the ground state energy cal- 
culated for a planer orbital (T points to the direction of 
d x 2_ y 2), which saddle point corresponds to the canting 
angle. 

Fig. H The bonding/anti-bonding splitting of the band 
by the inter-layer hopping. The small (large) splitting 
case is depicted as (a) ((b)). 

Fig. ^. Canting saddle point £ as a function of J5 (for 
the case £ c < 1). 

Fig. |t]. Schematic sketch of the calculated phase diagram 
and the correspondence with the experiments. 
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